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Abstract. The collision of a quantum Gaussian wave packet with a square barrier 
is solved explicitly in terms of known functions. The obtained formula is suitable for 
performing fast calculations or asymptotic analysis. It also provides physical insight 
since the description of different regimes and collision phenomena typically requires 
only some of the terms. 
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1. Introduction 

A paradigmatic textbook example pQ of a quantum scattering process is the collision of 
a particle of mass m, initially represented by a minimum-uncertainty-product Gaussian 
wave packet, with a "square barrier" potential, 



In spite of its simplicity, this model allows to observe and study a number of interesting 
quantum phenomena, such as tunneling OEllHElEllZllSlini, resonances [TU] . 
incidence- reflection and incidence-transmission interferences jTJ, the Hartman effect 
[121 HI IIHI UH EES! EE], and time delays or reflection due to a well when the height of the 
potencial, Vo, is negative [17] . It is also used as a standard model to check the validity 
of approximate numerical methods fElGH], or to exemplify and test different theories 
for temporal quantities such as arrival, dwell, or decay times, the asymptotic behaviour 
at long times [201 |2U 122] , and quantum transition state theories [23] . 

Except for approximate analytical treatments |12j . the solution of this time- 
dependent scattering model has been always obtained by numerical methods, typically 
using "Fast Fourier Transforms" , tridiagonal systems, or by linear combination of the 
solutions of the stationary Schrodinger equation. Goldberg, Schey, and Schwartz were 
the first to solve the model by means of a difference equation ^7j. While for many 
purposes these methods may be sufficient, for other applications the numerical approach 
may be cumbersome, time consuming or even useless. Therefore, an analytical solution 
in terms of known functions is of much practical interest. In general an analytical 
solution is not only useful as a fast computational tool; it is also valuable because of the 
physical insight than can be gained from it, very often by means of approximations that 
extract the dominant contributions in different limits, parameter ranges and regimes. 

In this paper we present an explicit expression in the momentum representation 
of the time dependence of a Gaussian packet incident on a square barrier in terms of 
known functions. There are a few other time dependent solvable scattering models: a 
Lorentzian state impinging on a delta function potential [2U 123 EE] or on a separable 
exponential potential [2E1 EZ] ; and a cutoff plane wave impinging on a square barrier 
[2EH2H]; or on a step barrier [30]. While all have been useful and illustrative of several 
scattering properties, the present model is the only one that combines simultaneously 
a local potential, resonances, and a physical Hilbert-space initial state with finite 
moments. 

2. Explicit solution 




otherwise. 



if - d/2 < x < d/2 



(1) 



The initial wave function is assumed to be a minimum-uncertainty-product Gaussian 
wave packet located in the left half-line and with negligible overlap with the potential. 
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In momentum representation it is given by 

(jpf\ij}(t = 0)) = ( —4 J exp \-S x (p - p c f /Ti 2 - ip'x c /h\ , (2) 

where x c and p c are the initial mean position and momentum respectively, and 5 X is 
the variance of the state in coordinate representation. Its time evolution can then be 
written in terms of eigenstates of H, \4> p i) jSDE21, as 

^(P,t) = / ( P \<f> P ')e- lE ' t/h (p'mt = 0)> dp', (3) 



with E' = p' 2 /2m. These states are given explicitly in coordinate representation by 

f l e ip ' x/h + Re- ip ' x/h x<-d/2 
{x\4>j,) = I Ce ip " x/h + De~ ip " x/n -d/2 < x < d/2 (4) 

[ Te ip ' x/h x > d/2, 

where p" = \Jp' 2 — 2mV§. The value of the coefficient / will be taken as 1. Note 
that these (delta-normalized) eigenstates are associated with an incident plane wave of 
momentum p' for p' > and with an outgoing plane wave with momentum p' for p' < 0. 
The coefficients R, C, D and T are determined by continuity of the wave function and 
its first derivative, 

sm(p"d/h)e- ip ' d / n 



R(p') 

C(p') 
D(p') 



2tt (p 1 ) 

(l _|_ £\ e -i(p'+ P ")d/2h 



V p" 



2VL (p') 
e~ 



2tt (p 1 ) 

p —ip'd/h 



where 



p" p' 



n (p') = cos(p"d/h) -l^ + ^j sm(p"d/h). (6) 

Fourier transformation of (x\(f) p >) in Q gives five terms proportional to the coefficients 
/, R, T, C, D respectively (generically denoted by A hereafter). Substituting these terms 
into (JHJ) with the initial state in Q, we obtain 



V<p, t) = irhh- 1 ' 2 / [ 9I (p') + g R (p') + g c (p') + g D (p) + g T (p)] e^dp' 

where r = {2-nTiy 1 ^ 2 (Ht) ^ • The term in the exponential <p(p') is given by 

-ip' 2 t 5 x (p'-p c ) 2 ip'x c ip'd 

Hp) = ^K 1? r~lh ] (8) 



Gaussian wave packet impinging on a square barrier 



4 



and 



9i(p) 
9n{p') 
9c V) 
gnip') 
9t(p') 



p — p' + iho + 
R(p')e ip ' d / h Q i P d / 2h 

p + p' + ihO + 
-2iC(p')e i P' d / 2h 

p — p" 
-2iD(p')e i P' d / 2h 

p + p" 
—T(p')e ip ' d ^ e _i P d / 2?i 



sm[{p-p")d/2h] 
sin[(p + p")d/2h] 



(9) 



p — p' — ihO + 

Note the three explicit ("structural" poles at p\ = p + ihO + , p' T = p — ihO + , 

and p' R = —p — ihO + , in addition to the poles of the functions C, D, R and T, p'j 
(j = 1, oo), which are zeros of fi(p')- All the poles lie in the lower half complex plane 
except for the incidence pole, p\. 

The integral in Eq. ((7j) can be solved by completing the square in Eq.(JSJ) and 
introducing the variable u as 
p' — z 



u = 



f 



with 



/_ [h 2+t 2mh 



-1/2 



u is zero at the saddle point z, defined by 



m 



Amp c 5l — (x c + d/2) Ti 2 t — i2mh [mS x (x c + d/2) + p c S x t] 



(10) 



(12) 



4m 2 6 2 + t 2 h 2 

and becomes real along the steepest descent path, a straight line with slope —th/(2mS x ). 
We now deform the contour and integrate along the steepest descent path, namely along 
the real u axis, 

V(p,t) = ifrhh- 1 ' 2 e-M/^W 



x / [9i (u) + 9r (u) + g c (u) + g D (u) + g T (u)] e u du, 



where 



V = 



' 2p c 5 x _ . {x c + d/2Y 

h 2 1 



h 



4 I + 



-1/2 



(13) 



(14) 



h z ' 2mh / 

gAiu) = gAip' — fu + -2), and Y u goes from —00 to +00 including a circle around the 
poles that have been crossed by the contour deformation. Since the gi, g^, gx and 
gcD = 9c + 9d are meromorphic functions with simple poles, it is useful to extract 
explicitly the singularities and leave the remainder as an entire function, h, 

Kr 
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9r(u) = + £ + h R {u) 

U-U R jr[U~ Uj 

9t(u) = + E + M«)- (15) 



Again the poles in the u complex plane may be separated into "structural" 
Uj = + i0 + - z) 

u R = r i (- P -io + -z) 

u T = f- 1 (p-iO + -z), (16) 
and "resonance" poles (see the Appendix), 

u j = f~ 1 (Pj-z) j = l,...,oo. (17) 

The residues TZ are given by 

Hi = - f- 1 e ipd/2h 
U R = f- 1 R{p' R )e ip 'n d/n e ipd/2h 
R^)e^i i/h e ipd l 2h F 3 



K 



Rj 



Cj 



n 



Dj 



[p + p'j + iho+) f 

-2iC(p' j )e ip ^ d/2n sm[(p - $)d/2h]Fj 
(p - p'^ f 

-2iD{p' j )e ip ^ d/2h sin[(p + p'^d/lh^ 



(p + p-)f 
H T = f~ l T(p' T )e ip T d/h e - ipd/2h 
T{p' 3 )e ip 'i dln e- ipd l 2h Fj 
(p ~ Pj ~ i%0 + ) f 



K 



Tj 



where 



dQ{p')/dp' p'=p'/ ^ 
Taking into account the integral expression of the w-function, w(z) = e~ z2 erfc(-iz) 



-u 2 



w( z ) = — f du, (20) 

27T Jr_ u — z 

where T_ goes from — oo to oo passing below the pole, and the relation w(—z) = 
2e~ z2 — w(z), the wave function may finally be written as 

V>(p,t) = ifrhh- 1/2 e-( 5xp2jh2 ) W {iirHMuj) - mTl R w{-u R ) 

oo 

- mTZ T w(-u T ) - in E \J^Rj + T^Cj + T^Dj + T^Tj]w(-Uj) 
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Figure 1. \ip(p, t)\ 2 as a function of momentum at t = 5 a.u., when the collision has 
been completed. The solid line corresponds to the exact solution, whereas the dots 
correspond to the reflection and transmission (structural) terms and three resonance 
pole terms from QH - Atomic units are used in all the calculations, and m = 1. p c = 20, 
d = 2.5, Vo = 188, S x = 100 and x c is located 50 atomic units to the left of the center 
of the barrier potential. 



which is the main result of this paper. The remaining Gaussian integrals are quite 
generally only a minor correction and may be evaluated with a few terms of the series 



which gives 

/ ^ Hu)e~ u du = ^ \h(u = 0) + £ ^ = 0) 

Note the form of the solution in (|2ip. There are structural terms associated with 
incidence, transmission and reflection (/, T, and R terms for short), resonance terms, 
and the /i-corrections of (j23J). 

3. Examples 

In many cases just a few terms in (|2T|) provide an accurate approximation to the exact 
result. If only the asymptotic behaviour, after the collision has been completed, is of 
interest, the structural terms R and T, associated with reflection and transmission, 
plus a few resonance pole contributions (again for reflection and transmission) are quite 
sufficient. The number of resonance poles that have to be added depends on the energy 
of the collision. In figure 1, for example, we show a case where only three resonance 
poles have been used in the R and T terms, and no correction term from Eq. has 
been included. 

There are also cases where only two terms are required during the collision: the 
collision with a very opaque barrier (R and I terms), see Fig. 2, and collisions for 



[h R (u) + h c {u) + h D (u) + h T {u)}e- u du\, 



(21) 




n=0 



(22) 
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Figure 2. \ip(p, t)\ 2 as a function of momentum, at t = 2.5 a.u., during the collision 
process. The solid line corresponds to the exact solution whereas the dots correspond 
to incidence and reflection structural terms of (|21[) only, and no resonance poles taken 
into account. p c = 20, d = 3, Vo = 400, 8 X — 100, and the initial wave packet is 
centered 50 atomic units to the left of the barrier center. 
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Figure 3. t)\ 2 as a function of momentum, at t = 2.7 a.u., during the collision 

process. The solid line corresponds to the exact solution, whereas the dots correspond 
to taking from 1|21|) incidence, reflection, transmission, and barrier (C and D) terms, 
plus three resonance pole terms. When one correction term is added, the 

approximate result is indistinguishable from the exact one. Same barrier and wave 
packet parameters as in Fig. 1 

very broad an energetic wave packets (J and T terms). This later case has been used 
recently to explain a transient, classically forbidden enhancement of the high momentum 
components plo"! EIo] as an incidence-transmission interference effect |TT] . 

In Figure 3, the case of a very general (typical) situation is depicted. Incidence, 
reflection, transmission, and barrier (C and D) terms from (j21J) . plus three resonance 
pole terms provide a fairly good approximation for the wave function at an intermediate 
time during the collision, although the approximate solution clearly overestimates the 
exact value at the reflection and transmission peaks, and for values of momentum 
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between 20.3 and 21 a.u.. The inclusion of the first correction term in (J23)) leads to 
a curve which is indistinguishable in the scale of the figure from the exact result. 

Let us also note that pijl is suitable for asymptotic analysis (e.g. short and large 
times) by means of the asymptotic series of the w-function, see e.g. (23 HZ] • The w- 
functions may be considered the elementary propagators of the Schrddinger transient 
modes [37] and play a prominent role in the explicit solution given in (}2T]) . 
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Appendix A. Calculation of the resonance poles 

The resonance poles are obtained by solving the transcendental equation Q(p') = 0. A 
systematic way to do it is based on the differential equation [3*H] 

dp' _ idV Q y/m - p' ^ 



dV 14 (£-2)- 
As boundary condition for the integration of (jA.l|) we take p' = at 



2 2 

n 7T 



Vo = ~^ ] n = 0,±l,±2,... (A.2) 

The integration of (jA.l)) has to be carried out until the value of Vq corresponds to the 
physical barrier we want to examine. The integration however cannot be performed 
along the real axis because of the factor Vq in the denominator of (jA.lj) . To avoid 
the singularity at Vq = 0, an imaginary term is added to the potential, analytically 
continuing the differential equation. A parabolic path has been used across the complex 
Vq plane from the initial values of Vq in ()A.2j) to the final real value of Vq. 
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